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Abstract 



We have computed E>k with staggered fermions, using two different methods. The 
numerical simulations were performed on a 16 3 x 40 lattice in full QCD with @ = 5.7. 
' We also tried an improved wall source method in order to select only the pseudo- 

Goldstone bosons and compare the numerical results obtained with those from the 
conventional wall source method. We have studied a series of non-degenerate quark 
anti-quark pairs and saw no effect on Bk, although effects were seen on the individual 
■ terms making up Bk- 

ON 

1 Introduction 

In the standard model, there are two kinds of CP violation: the indirect CP violation and 
the direct CP violation. [[|, |j, H|. The indirect CP violation comes from the fact that 
electro-weak K°-K° mixing causes the neutral kaon eigenstates not to respect CP symmetry 
||. The direct CP violation comes from the CP violating effective operators (e.g. penguin 
operators) @, ||, which can directly contibute to the K decay amplitude. The indirect 
(direct) CP violation is parametrized by the phenomenological quantity e, while the direct 
CP violation is parametrized by e' [p], 0, |3], |J . Both e and e' can be measured experimentally. 

The low energy effective Hamiltonian of the electro- weak interaction is derived by de- 
coupling heavy particles such as W and Z bosons and the t, b, c quarks in the stan- 
dard model H, [7|, gfl. The coefficients in this effective Hamiltonian are functions of 
the Cabibbo-Kobayashi-Maskawa flavor mixing matrix elements Vqkm- This low energy 
effective Hamiltonian includes a AS = 2 term which belongs to the (27,1) representa- 
tion of the SU(3)l <8> SU(3)r flavor symmetry group and to the (1,0) representation of 
the SU{2) L g) SU(2) R isospin symmetry subgroup. This AS = 2 electro-weak effective 
Hamiltonian causes the K°-K° mixing in the standard model |3|, |6|, 0, ||. Therefore, the 
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AS = 2 effective Hamiltonian is connected to the e parameter. Using this technique, e can 
be expressed in terms of the standard model parameters as follows: 



e = CB K F(Vckm; 
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B K = Z(a a (ji))B K <ji) . 



(2) 

(3) 
(4) 



Here Bk is defined as a renormalization-group invariant quantity, and the function 
F(Vckm, ™>t/Mw, • • ■) is given in Ref. Jl|, |], £§. Once Bk is determined theoretically, 
we can narrow the domain of | Vtd | and the top quark mass |l|], using the experimental 
determination of e. 

There have been a variety of methods used to calculate Bk, including chiral perturba- 
tion theory, hadronic sum rules, QCD sum rules, l/N c expansion and lattice gauge theory. 
Lattice gauge theory has the virtue of making the smallest number of assumptions and is 
exactly equivalent to QCD in the limit of infinite volume and vanishing lattice spacing. For 
this reason, we have adopted lattice gauge theory to obtain B K . 

In order to calculate B K in lattice gauge theory, one needs to find a set of operators which 
can describe on the lattice the same physics as the continuum AS = 2 four-fermion operator. 
There have been two methods to implement fermions in lattice gauge theory: the staggered 
fermion method and Wilson fermion method. For the weak matrix elements involving the 
pseudo-Goldstone bosons, it is very useful to take advantage of the exact Ua(1) symmetry 
of the staggered fermion action, which is not manifest in Wilson fermion action ||, 1Q| . We 



have used the staggered fermion method to calculate B K because of this advantage. 

There are two methods to transcribe the continuum AS = 2 weak matrix elements to the 
lattice with staggered fermions R II]: the one spin trace formalism and the two spin trace 
formalism. The four fermion operators can be expressed as products of operators bilinear in 
the fermion fields. In the one spin trace formalism, each external hadron is contracted with 
both bilinears of the four fermion operator simultaneously. In the two spin trace formalism, 
each external hadron is contracted with only one of the bilinears in the four fermion operator. 
Until now, the two spin trace formalism (2TR) has been used exclusively in calculations of 

and jr 
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weak matrix elements with staggered fermions, |T2|, [L3], [L4J], |T5], [16], [T7|, [T 
Recently, the one spin trace formalism has been developed to a level which permits it to be 
used for numerical simulations of weak matrix elements such as Bk JTTI] . We have used both 



formalisms to calculate Bk and the results are compared in this paper. 

Lattice calculations introduce their own systematic artifacts and errors which can be 
sizable. One of the principal sources of systematic errors is finite volume. The results of 
a finite volume comparison were reported in Ref. fL6f, where it was argued that the finite 



volume effects on Bk are quite small when V x T > 16 x 40 at (3 = 6.0 (a ~ 2 GeV), 
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Nf = 0. Another source of systematic errors comes from the quenched approximation, 



neglecting all of the internal quark loops. Quenched effects were studied in Ref. fllSfl , which 



concluded that the effect of quenching can not be large. Another systematic error which is 



called scaling violation comes from the finite lattice spacing (a ^ 0). In Ref. fUfl , the scaling 
corrections were argued to be of a 2 order. In other words, 0(a) corrections do not exist at 
all for the staggered fermion operators of Bk- Related to possible scaling violations is a 
possible discrepancy between gauge-dependent Landau gauge operators and gauge invariant 
operators. The Landau gauge operator implies that the quark anti-quark propagators are 
fixed to Landau gauge and that gauge links between the quark and anti-quark fields are 
omitted. The question was raised as to whether the Landau gauge operators might cause 



the large scaling violation which had been noticed originally in Ref. (T^|. However, it was 
reported in Ref. jn| that the results of both Landau gauge operators and the gauge invariant 
operators were numerically found to be consistent with each other. 

The purpose of this paper is to report and analyze the numerical data for Bk as well as its 
individual terms and to interpret the numerical results in terms of various physical models. 



Part of the preliminary results have already appeared in Ref. [21]. In this paper we will 
address the following five issues through the interpretation of our numerical results. The first 
issue is how to select the pseudo Goldstone boson state exclusively. For hadron spectrum 
measurements, the sink operator picks up a specific hadronic state exclusively. In contrast 
to hadron spectrum measurements, the operators for weak matrix element measurements 
do not select a particular hadronic state. We need to impose a symmetry requirement on 
the wall source such that all the unwanted states are decoupled in the weak matrix element 
measurement. We have tried an improved wall source method (called cubic wall source f22|| ) 
to do this. The second issue is whether the numerical results in the one spin trace form are 
in agreement with those in the two spin trace form. Theoretically the difference between the 
two formalisms is supposed to vanish in the limit of a —>■ |L1] for Bk- We have tried both 
the one spin trace formalism and the two spin trace formalism to calculate Bk- The results 
are compared in this paper. The third issue is the effect of non-degenerate quark antiquark 
pairs on Bk and the individual components making up Bk- The kaon is composed of s and d 
valence quarks. Here a non-degenerate quark anti-quark pair implies that the s valence quark 
mass is different from the d valence quark mass, whereas a degenerate quark anti-quark pair 
implies that both valence flavors have the same mass. The effect of non-degenerate valence 



quarks on Bk in quenched QCD was mentioned briefly in Ref. ]13| where a small difference 
of only marginal significance was found. In this paper, we investigate in detail the effects of 
non-degenerate valence quark anti-quark pairs on Bk and its individual components. From 
the theoretical point of view, one effect of non-degenerate valence quark anti-quark pairs 
can be related to the rf hairpin diagram in (partially) quenched QCD |23| , |24| , |25| . The 
forth issue is whether quenched chiral perturbation theory is compatible with the numerical 
results of Bk |[23| . Quenched chiral perturbation theory also predicts the chiral behavior of 
the individual terms making up Bk [fffil - It is good to know how reliable these theoretical 



predictions are numerically. The final issue is whether there is any dynamical fermion effect 



on Bk- This question was addressed originally in Ref. ||15|| . We will re- visit this question and 
see how important the internal fermion loops are to Bk- It is important to understand the 
difference between full QCD and quenched QCD in Bk both theoretically and numerically. 
This paper is organized as follows. In section 2, we will describe the technical details 
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in doubling the lattice for quark propagators and explain the improved (cubic) wall source 
method in a self-contained manner. In Section 3, we will specify the lattice operators for 
Bk in brief while leaving the details to adequate references. In Section 4, we describe the 
parameters for the gauge configurations we generated and the measurement parameters for 
Bk- In Section 5, we present the numerical data for M.r (numerator of Bk) and the vacuum 
saturation amplitude (denominator of Bk)- The main emphasis is put on the determination 
of the plateau region. We also discuss the improved wall source with the wrong flavor 
channel, which is supposed to vanish in the limit of a — > 0. In Section 6, the numerical 
results from the data analysis are interpreted. The improved (cubic) wall source method is 
compared with the conventional method. The data of the one spin trace form and of the two 
spin trace form are compared. Enhanced chiral logarithms in the individual terms making 
up Bk appears to be seen numerically as well as the additional divergence which arises from 
the effect of the non-degenerate quark anti-quark pairs. We emphasize the large effect of 
non-degenerate quark anti-quark pairs on the individual components of Bk- We compare 
our data with other groups (Kilcup, Sharpe and Ukawa et al). We describe the covariance 
fitting procedure for Bk and the negligible effect of non-degenerate quark anti-quark pair 
on Bk- Section 7 contains a brief summary and our conclusions. 



2 Quark Propagators and Wall Sources 

Here, we explain the technical details of the quark propagators and the improved wall source 
method. 

Because the lattice size in the time direction is finite, the pions propagating around 
lattice in the time direction can in principle contaminate the measurements of Bk- There 
have been two proposals for avoiding this contamination. The first proposal is to impose 
Dirichlet boundary conditions in time on the quark propagator and to place the wall source 
near the boundary |H| [TJ, [H| . The second proposal is to double the lattice along the time 



direction and to use periodic boundary conditions in time for the quark propagator [10], [13 
Dirichlet boundary conditions cause a certain number of time slices near the wall source to be 
contaminated by reflections off the boundary. The time slices lost due to reflections overlap 
with those contaminated by p mesons. Lattice doubling in the time direction suppresses 
the backward propagating pions which must travel over the whole lattice size in the time 
direction before they can contribute to the measurements. In contrast to Dirichlet boundary 
conditions, lattice doubling does not introduce any reflection from the boundary. Hence, the 
signal from lattice doubling is much cleaner than that obtained using Dirichlet boundary 
conditions. Unfortunately, lattice doubling takes twice the computational time since the 
lattice doubling needs two undoubled quark propagators for a given wall source (forward 
and backward) instead of the one quark propagator required when using Dirichlet boundary 
conditions. 

We adopt the second solution of doubling the lattice for our numerical simulation of 
Bk- Here, we describe this second method in detail. The periodic and anti-periodic quark 
propagators in time on the undoubled lattice are: 

52(D + M) {gtX yGp(x,y) = 6 ZiV (5) 



4 



and 

J2(D + M) (z , x yG A (x,y) = 5 z , y , (6) 

X 

where (D+M) represents the Euclidean Dirac operator, and Gp(x, y) and G A (x, y) represent 
the Green's functions with periodic and anti-periodic boundary conditions on the undoubled 
lattice. The source h(y) is introduced in the numerical simulation as follows: for < x t < N t , 

X(x) = ~ J2(Gp(x, y) + G A {x, y)) ■ h(y) , (7) 
A y 

X (x + Nj) = \ y) - G A {x, y)) ■ h{y) (8) 

1 y 

where t is the unit vector in time direction and < yt < N t (N t is the undoubled lattice size 
in the time direction). The xi x ) field in the above satisfies the following: for < Xt < N t , 

J2(D + M) {z , xyX (x) = h(z) (9) 

X 

and for N t < x t < 2N t , 

J2(D + M) M ■ xix + Nj) = . (10) 

X 

We use periodic and anti-periodic boundary conditions on the undoubled lattice and take the 
sum and the difference to obtain a quark propagator which is periodic on the doubled lattice 
in the time direction. Note that the fermionic field xi x ) is periodic on the doubled lattice 
in the time direction. One might ask whether periodic boundary conditions in the time 
direction make any difference compared to the anti-periodic boundary conditions, which are 
the physical ones at a finite temperature. The answer is that periodic boundary conditions 
make no difference as long as the volume is large enough. The physical eigenstates in the 
confined phase of QCD are hadrons, not quarks nor gluons. In the confined phase, a quark 
must be confined with other quarks or anti-quarks within a small volume of typical hadronic 
size (~ 0(lfm)) before it can acknowledge the existence of the boundary with the volume 
much larger than the hadronic size. 

There have been a number of attempts to enhance the overlap with the lightest particle 
compared to that of the excited states so that one can see the asymptotic signal (exp(— M | 
t I)) at smaller time separations and over the longer plateau |TH EEJ. The wall source 



method has an advantage because it enhances the signal of the hadron propagators with 
respect to the point source method. For hadron mass spectrum measurements, the sink 
operator possesses the symmetry of a specific hadronic state. In contrast to hadron spectrum 
measurements, since the operator in the electro-weak effective Hamiltonian does not select 
any particular hadronic state by itself, the weak matrix element measurements require the 
symmetry properties of the wall source to determine the specific hadronic state. There have 
been two attempts to improve the wall source such that it can exclude contamination from 



unwanted hadronic states: the even-odd wall source method |26| and the cubic wall 



source method originally proposed by M. Fukugita, et al. in Ref. 
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Here, we explain the cubic wall source method in detail. The Bx measurements requires 
that the pseudo-Goldstone Kaon mode should be selected exclusively. Hence, we will restrict 
our discussion to the cubic wall source operator of pseudo-Goldstone mode. Let us start with 
definitions and notation. The symbol a represents one of the 8 vertices in the unit spatial 
cube {i.e. a 6 {0, l} 3 ): 

a e {(0,0,0), (1,0,0), (0,1,0),-. -,(1,1,1)} . (11) 

We define W a {n) as follows: 

W- a (fi) = <k2m+a • (12) 

meZ 3 

The cubic wall source operator for the pseudo-Goldstone pions can then be expressed as 
follows 

^source = E x'fa * = 0)W a {n) e(a)5 , 6 Wt{n') X \n\ t = 0) , (13) 

n,n',a 

where e{x) = (— \Y 1+X2+X3+Xi , a anc j ^ are co i or indices and \ is a staggered fermion field. As 
an example, let us choose the sink operator to be the bilinear operator with spin structure 
S and flavor structure F: 

O sink (y) = x{yA){ls®tF) AB x{yB) ■ (14) 

Ua = 2y + A, where y G Z 4 is the coarse lattice coordinate and A e {0, l} 4 is the hypercubic 
coordinate. Then the correlation function is 

\^sink(?y) ^source) 

= E (75 ® Z F ) AB G c > a (y B , n)W a (n)e(z)5 a , b W a (n')G b > c (n', y A ) 

n,n',a 

= E (7575 ® MF) AB G c ' a {VB, n)W a {n)5 a> M{n')[G c ' b {y A , ft')}* 

n,n',a 

= (7575 i Mf) ab E i>i(yB)bPi{yA))* (is) 

c,6,a 

where G{x, y) = \[Gp{x, y) + G A {x, y)] for the domain < x t , y% < 2N t and 

^(y A ) = Y.G c > b (y A ,n)W- a (n) ■ (16) 

n 

Here ip a b {y A ) is actually what we calculate on the computer using a conjugate gradient 
method with the wall source set to W a . For the fx and the vacuum saturation part of 
Bx, Eq. ( Jl5| ) is used in our numerical simulation. The idea of the above example (bilinear 
sink operators) can be extended as a whole to the four fermion operator measurements (for 
example, Bx) without loss of generality. 
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3 Operators Computed 



Here, we present a set of lattice operators which describe the same physics as the continuum 
AS = 2 operator and which have been used for our numerical calculation of B&. In the 
continuum, B K is defined as 



B M = 8/ r >m ... .. jinWnl =1 -. -7 1 T^m ( 17 ) 



(if 


S7m(! _ 75)^ «7m( 1 _ 75)^ 
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s7m75^ 
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K°) . 



The numerator of is a matrix element of the AS = 2 four-fermion operator with the 
neutral K meson (kaon) states. The denominator of Bk represents the vacuum saturation 
approximation of the numerator, which inserts the vacuum state between the two bilinears 
of the AS* = 2 four-fermion operators. For the denominator, one needs to prescribe a lattice 
bilinear operator which corresponds to the continuum axial current: 



A-v = Xs(7m5 ® &)Xd ■ (1* 



Here the notations are the same as those adopted by Ref. JjJ [10], [TTH . The axial current is 
chosen such that the flavor structure is identical to that of the pseudo Goldstone kaon of the 
exact Ua(1) symmetry in the staggered fermion action. 

For the numerator, we need some particular set of lattice four-fermion operators which 
correspond to the continuum AS" = 2 operator. There are two methods to transcribe the 
continuum AS = 2 operator on the lattice with staggered fermions: the one spin trace 
formalism and the two spin trace formalism |Tt], [TJJ. In the continuum, there is no 
difference between the one spin trace form and the two spin trace form of the AS = 2 
operator, since they are connected by the Fierz transformation. However, in the staggered 
fermion method, the one spin trace form is different from the two spin trace form due to a 
pure lattice artifact of 4 degenerate flavors. 

In the two spin trace formalism, the four-fermion operator in the numerator in Eq. ( [P7D 
is transcribed to the lattice || |10[ as a sum of four operators: 

Offit = (VxP)%&+(VxP)%* 

+ (A x Pf™ a + (A x P)2« (19) 

where V (or A) represents the vector (or axial) spin structure, P represents the pseudoscalar- 
like flavor structure and the subscript ab; ba (or aa; bb) represents the color indices of the 
quark fields (the details of these notations are described in Ref. ||11||). The operators in Eq. 
(|i~9l) have the same chiral behavior in the limit of vanishing quark mass as the continuum 
AS = 2 operator JTD|, [TTJ, ^J. In addition, the matrix elements of O^r show the same 
leading logarithmic dependence on the renormalization scale as the continuum AS = 2 
operator [TT|, |3T | . 



In the one spin trace formalism, the four-fermion operator of the numerator in Eq. (|T 
is transcribed to the lattice as follows: 

0™ = (^p)™ + (yxP)3 

+(AxP)g + (AxPr 



1 aa\b 

+ 1TR (20) 

T ^chiral partner " V w / 
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where again the details of this notation is given in Ref. [11]. In contrast with the two spin 
trace formalism, the individual terms in Eq. (f20| ) do not possess the same chiral behavior 
as the continuum AS = 2 operator |IT| . We must add Cchimi partner m or der to preserve the 



correct continuum chiral behavior pl| . By imposing the correct chiral behavior on O^r 



we determine the chiral partner operator JTT| as follows 



,1TR 
)aa;bb 



^chiral partner _ (V" X S)^.^ a + (V X S)] 

+(A x S)l™ + (A x S)ZS b ■ (21) 

This forces the resulting operator to respect the continuum chiral behavior. The next ques- 
tion is whether the whole operators including the additional chiral partner operators still 
have the continuum leading logarithmic behavior. First of all, we need to choose the basis 
operators such that they belong to the identity representation with respect to the 90° axial 
rotation group (a subgroup of the exact Ua(1) symmetry group). This particular choice of 
the basis operators guarantees the analytic chiral behavior of the continuum. Second, we 
find that an eigen-operator (Eq. (|20| ) and Eq. (|2T|) ) possesses the same leading logarithmic 
behavior as the continuum AS = 2 operator [[□]]. For Bk measurements, the difference 
between the one spin trace operators and the two spin trace operators vanishes as a — ► 



11]. We have used both one spin trace and two spin trace operators to calculate B K on the 



lattice and the numerical results are compared later in this paper. 



4 Simulations and Measurement Parameters 

In this section we describe the simulation parameters for the gauge configurations and Bk 



measurements. Our old hadron mass calculation [27 on a 16 3 x 32 lattice at (3 = 5.7 informed 
us that a longer time dimension allows a more precise fitting of the hadron propagators' 
exponential time dependence. For this reason the volume of the configuration was chosen 
as 16 3 x 40. The coupling constant was (3 = 5.7 (1/a m 2 GeV). The dynamical effects of 
two degenerate flavors of staggered fermions with a mass 0.01 were incorporated into the 
gauge configurations, using the hybrid molecular dynamics R-algorithm |28| . The sea quark 
mass (0.01) has been fixed through all the measurements even though various valence quark 
masses were chosen for the Bk measurement. The gauge configurations were updated by the 
hybrid molecular dynamics R-algorithm with molecular dynamics step size 0.0078125 p8 
and a trajectory length of 0.5 time units. 

Now the measurement parameters for Bk'- Every 60 trajectories Bk has been measured. 
The total number of the gauge configuration samples for Bk measurements was 155. We have 
used both cubic wall source and conventional even-odd wall source methods to create the 
pseudo Goldstone boson. We used two separate wall sources to create K° and K° mesons. 
The distance between these two separate wall sources was 36 lattice units. For each Bk 
measurement, both wall sources were shifted by 15 lattice units in the time direction from 
the position used in the previous measurement, while the distance between them was fixed 
to 36 through all the measurements. The valence (quark, antiquark) mass pairs for the K 
mesons were (0.01, 0.01), (0.02, 0.02), (0.03, 0.03), (0.004, 0.01), (0.004, 0.02), (0.01, 0.03) 
and (0.004, 0.05). The quark propagators were calculated, using the lattice doubling method 
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described in the previous section. The stopping condition of the conjugate gradient residual 
for the quark propagator was set to 1.0 x 10 -8 . The quark propagators were gauge fixed to 
Landau condition with a numerical precision of a /L g 2 (d fJi A fJi ) 2 < 1.0 x 10 -7 . 

5 Data Analysis 

We now present the numerical results of Mk (the numerator of Bk) and the vacuum satu- 
ration amplitude Ml (the denominator of Bk) with respect to the lattice Euclidean time 
for Kaon with a quark anti-quark mass pair = (0.01, 0.01). Meanwhile, we explain how the 
central plateau region has been selected to determine Bk- We also discuss the numerical 
results for unrenormalized (naive, bare) Bk with respect to the various quark and anti-quark 
masses. In addition, we will present the numerical measurement of the wrong flavor channel 
((V + A) (gi S) 2TR in order to see how much contamination comes from the operator mixing 
and the excited hadronic states which are supposed to vanish in the limit of a — > 0. 

5.1 The Denominator: Vacuum Saturation Amplitude 

We define the denominator of Bk as the vacuum saturation amplitude: 

Ml = \{K G \ Sl,l,d | 0)(0 | s ltil5 d | K°) . (22) 

Numerical data for Ml is presented in Figure |I| (the even-odd wall source method) and in 
Figure (the cubic wall source method). The data points in Figure [I], 0, and |3] are obtained 
by single elimination jack-knife method. Each data point in Figure |l|(the even-odd wall 
source) has about twice larger error than that in Figure |2] (the cubic wall source) . However, 
one needs to notice that the cubic wall source method takes four times longer time to compute 
than the even-odd wall source method. In Table |], the results of the covariance fitting of the 
vacuum saturation together with their x 2 per degree of freedom are given with respect to the 
various fitting ranges in the lattice Euclidean time for quark mass m q a = 0.01. From Table 
[L], for the cubic wall source method the minimum of the x 2 per degree of freedom extends to 
the fitting range of 11 < t < 24. For the even-odd wall source method, the minimum of the 
X 2 per degree of freedom occurs in the fitting range of 13 < t < 23. It is important to keep 
in mind the fact that both the one spin trace form and the two spin trace form of Bk have 
the vacuum saturation in common. 

5.2 The Numerator: Mr 

We define the numerator of Bk as: 

M K = (K° | sj^l - l5 )d s Tm (1 - l5 )d | K°) (23) 

Numerical data for Mk calculated in the two spin trace formalism with the even-odd wall 
source method is presented in Figure §. Numerical data for Mk calculated in the two spin 
trace formalism with the cubic wall source method is shown in Figure Numerical data 
for Mk calculated in the one spin trace formalism with the even-odd wall source method 
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Fitting Range 


Cubic Wall Source 


Even-Odd Wall Source 


M V K 


X z /(dof) 


M V K 


X z /(dof) 


14 < t < 21 


12.17(66) 


0.40(63) 


49.2(30) 


0.40(68) 


13 < t < 22 


12.10(58) 


0.37(52) 


49.3(25) 


0.33(53) 


12 < t < 23 


12.15(60) 


0.37(48) 


48.7(24) 


0.61(65) 


11 < f < 24 


12.10(63) 


0.45(48) 


49.0(25) 


0.68(56) 


10 < t < 25 


11.76(59) 


1.06(51) 


47.0(27) 


1.24(76) 


9 < i < 26 


11.70(59) 


1.27(66) 


46.7(27) 


1.34(70) 


8 < i < 27 


11.61(58) 


1.31(62) 


45.2(24) 


1.40(68) 


7 < i < 28 


11.20(58) 


1.65(70) 


44.2(24) 


1.52(70) 


6 < t < 28 


10.74(60) 


1.72(65) 


42.9(24) 


1.47(67) 



Table 1: Here, we present the numerical results for the vacuum saturation amplitude (A4]Q 
with the quark mass pair (0.01, 0.01), calculated both in the even-odd wall source method 
and in the cubic wall source method. All the values in the table have been obtained through 
the covariance fitting to a constant on the bootstrap ensembles. 



is drawn in Figure |^. All the data points in Figures f|, and ^] are obtained by the single 
elimination jack-knife method. Each data point in Figure ^| (the even-odd wall source) has 
an about twice larger error than that in Figure [| (the cubic wall source). The nearby data 
points in Figure |3] and Figure [| (the even-odd wall source) have larger fluctuations than those 
in Figure § (the cubic wall source), which are reflected on x 2 /d.o.f. in Table |[ However, 
this was not observed in the results for vacuum saturation amplitude. It is important to 
note the fact that the computational time for the cubic wall source method is four times 
longer than that for the even-odd wall source method. In Table |2], the covariance matrix 
fitting results of M.k and its x 2 /(d.o.f.) are collected with respect to the fitting range in the 
lattice Euclidean time. From Table [2] it is difficult to choose the optimal fitting range as all 
of the x 2 /(d.o.f.) values are within 0.5 ~ 1 a for M.k calculated in the two spin trace form 
with the cubic wall source. This is also true for M.k calculated in the one spin trace form 
with the even-odd wall source. For these cases, we therefore choose the optimal fitting range 
consistent with the vacuum saturation amplitude AiJ^- 

For M.R calculated in the two spin trace form with the even-odd wall source, we notice 
that the optimal fitting range is 13 < t < 22. The M.k results in the cubic wall source 
method has more correlation (less fluctuation) between neighboring time slices than those in 
the even-odd wall source method, while this is not obvious for the measurements of vacuum 
saturation amplitude M^- 

5.3 The Ratio: B K 

In the previous sections, we have discussed the covariance matrix fitting result of M.r and 
Ai^ as a function of the fitting range. From this analysis of the fitting ranges, we determine 
the optimal fitting range for B^- The optimal fitting ranges we have chosen to determine 
Bk with respect to various quark mass pairs are summarized in Table [5| Once the fitting 
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Fitting 
Range 


Two Spin Trace Form 


One Spin Trace Form 


Cubic Wall Source 


Even-odd Wall Source 


Even-odd Wall Source 


M K 


xV(d.o.f.) 


M K 


X V(d.o.f.) 


M K 


x7(d.o.f.) 


14 < t < 21 


7.20(48) 


1.11(99) 


29.8(15) 


1.05(102) 


33.4(21) 


1.45(106) 


13 < t < 22 


7.23(29) 


0.88(87) 


29.9(15) 


0.87(83) 


33.6(17) 


1.18(82) 


12 < t < 2d 


7.26(28) 


0.75(75) 


28.7(14) 


1.47(80) 


34.0(17) 


1.19(68) 


11 < t < 24 


7.19(28) 


0.73(67) 


28.6(13) 


1.24(70) 


34.1(18) 


1.03(62) 


10 < t < 25 


7.16(29) 


0.66(60) 


28.1(14) 


1.29(59) 


33.7(16) 


0.96(58) 


9 < t < 26 


7.11(30) 


0.62(51) 


28.1(13) 


1.14(54) 


33.7(16) 


1.04(49) 


8 < t < 27 


7.11(30) 


0.60(52) 


27.7(12) 


1.27(61) 


33.6(16) 


0.95(47) 


7 < t < 28 


7.04(27) 


0.57(54) 


27.1(11) 


1.47(52) 


33.9(16) 


0.95(43) 


6 < t < 29 


6.97(26) 


0.61(50) 


27.3(12) 


1.52(54) 


33.8(15) 


0.90(45) 



Table 2: Here, we present the numerical results for the M.r with the quark mass pair (0.01, 
0.01), calculated both in the even-odd wall source method and in the cubic wall source 
method. All the values in the table has obtained through the covariance matrix fitting to a 
constant on the bootstrap ensembles. 



range is chosen, Bk can be determined through the covariance fitting to a constant. One of 
our fitting procedures is naive jack-knife analysis of the data (conventional) and the other 
uses the covariance fitting on the bootstrap ensembles. 

The numerical results for Bk calculated in the two spin trace form with the even-odd 
wall source method are shown in Figure |6|. The numerical results for Bk calculated in the 
two spin trace form with the cubic wall source method are drawn in Figure [7[ The numerical 
results for Bk calculated in the one spin trace form with the even-odd wall source method 
appear in Figure ||. 

Table |] is a collection of the numerical results of K meson mass and Bk with respect to 
the quark mass pairs calculated in the optimal fitting range by the jack-knife method. The 
results for K meson mass (tuk) are obtained by analyzing the results for axial current matrix 
element with an external K meson state (0 | ^(7075 ® £,§)D \ K), which are also used to 
obtain the vacuum saturation amplitude. We fitted the logarithm of the axial current data 
to the linear function of the Euclidean time A + m^t in the optimal fitting range. We would 
like to thank Pavlos Vranas for checking these mass results using his own fitting program. 
Table [| presents the results of Bk obtained by the covariance fitting over the optimal fitting 
range on the bootstrap ensembles. 

5.4 The Wrong Flavor Channel 

Here we address two important questions on the validity of our approach to Bk- The first 
question comes from the fact that the higher-loop radiative correction of the four-fermion 
composite operators cause the violation of the continuum spin and flavor symmetries. It 
is important to know non-perturbatively, how large is the contribution of such symmetry 
violating terms to the weak matrix element measurements for finite lattice spacing. Note, 
such terms are supposed to vanish in the limit of a — > 0. The smaller the contribution 
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Quark Mass 
Pair 


Fitting Range 


Two Spin Trace Form 


One Spin Trace Form 


Cubic Wall Source 


Even-odd Wall Source 


Even-odd Wall Source 


(0.004, 0.01) 


9 < t < 26 


13 < t < 22 


13 < t < 22 


(0.004, 0.02) 


11 < t < 24 


13 < t < 22 


13 < t < 22 


(0.004, 0.05) 


14 < t < 21 


15 < t < 20 


15 < t < 20 


(0.01, 0.01) 


11 < t < 24 


13 < t < 22 


13 < t < 22 


(0.01, 0.03) 


11 < t < 24 


13 < t < 22 


13 < i < 22 


(0.02, 0.02) 


13 < t < 22 


13 < t < 22 


13 < i < 22 


(0.03, 0.03) 


13 < t < 22 


13 < t < 22 


13 < t < 22 



Table 3: Here, we present the optimal fitting range with respect to quark mass pairs. The 
optimal fitting range implies that we can get the same average with smaller error bar for 

X 2 < i-o. 



Quark Mass 
Pair 


m K 


unrcnormalizcd Bk 


Two Spin Trace Form 


One Spin Trace Form 


Cubic Source 


Even-odd Source 


Even-odd Source 


(0.004, 0.01) 


0.219(2) 


0.557(32) 


0.547(54) 


0.656(73) 


(0.004, 0.02) 


0.277(2) 


0.641(32) 


0.630(39) 


0.698(58) 


(0.004, 0.05) 


0.406(2) 


0.731(34) 


0.687(44) 


0.759(88) 


(0.01, 0.01) 


0.253(2) 


0.600(27) 


0.579(33) 


0.688(46) 


(0.01, 0.03) 


0.348(2) 


0.710(26) 


0.689(30) 


0.748(37) 


(0.02, 0.02) 


0.347(1) 


0.709(25) 


0.689(28) 


0.746(33) 


(0.03, 0.03) 


0.421(1) 


0.768(23) 


0.753(26) 


0.781(28) 



Table 4: K meson mass and unrenormalized Bk versus quark mass pair: the Bk data are 
analyzed by the single-elimination jack-knife method over the optimal fitting range. The K 
meson mass, rriK was obtained by analyzing the results for (0 | 5'(7o75 <8> £,b)D \ K). 



Quark Mass 
Pair 


Unrenormalized Bk 


Two Spin Trace Form 


One Spin Trace Form 


Cubic Source 


Even-odd Source 


Even-odd Source 


(0.004, 0.01) 


0.546(23) 


0.548(42) 


0.675(63) 


(0.004, 0.02) 


0.627(28) 


0.663(37) 


0.700(47) 


(0.004, 0.05) 


0.732(30) 


0.715(40) 


0.771(78) 


(0.01, 0.01) 


0.595(22) 


0.607(34) 


0.682(45) 


(0.01, 0.03) 


0.717(23) 


0.708(31) 


0.771(37) 


(0.02, 0.02) 


0.727(24) 


0.707(29) 


0.773(36) 


(0.03, 0.03) 


0.791(21) 


0.758(26) 


0.804(25) 



Table 5: Unrenormalized Bk for each quark mass pair. The Bk data is analyzed by the 
covariance fitting over the optimal fitting range on the bootstrap ensembles. 
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from these wrong flavor channels is, the more reliable our connection between the continuum 
and lattice observables (either at tree level or at one loop level). The second question is 
how exclusively we can select the pseudo-Goldstone mode through our improved wall source 
method. In other words, we want to ask how reliably our wall source technique suppresses 
the unwanted hadronic states. 

We have chosen the ((V+A)(3S) 2TR operator in order to address the above two questions. 
The matrix element of this operator with external K mesons is supposed to vanish in the 
continuum limit (a — ► 0) of lattice QCD, due to the vanishing flavor trace. The numerical 
results for this wrong flavor channel are shown in Figure |9| for the even-odd wall source 
and Figure [10] for the cubic wall source. From Figures and [10, we notice that the value 
of the wrong flavor channel is extremely suppressed (less than 1% of Bk) in both cases. 
This implies that the unwanted operator mixing of ((V + A) <g> S*) 2 ™ should be at most 
1% of B K since it is suppressed by a s /(47r) as well as by the vanishing flavor trace. This 
fact that the unwanted operator mixing is smaller than 1% of Bk is of great significance 
to our approximate matching between the continuum and lattice composite operators. So 
far, one has neglected those terms of wrong flavor channels which enters at order g 2 , when 
one connects the lattice Bk to the continuum Bk at one loop level [K| [11, 19], [HJ. The 



main reason was that these wrong flavor channels will not contribute to Bk at all as a — > 0. 
Hence, the remaining difficulty was to know how large is the contribution of those wrong 
flavor channels at finite lattice spacing. Our non-perturbative measurements of one of wrong 
flavor channels shows that the contribution from these wrong flavor channels will be at most 
1% of Bk and so much smaller than the statical and other systematic errors at j3 = 5.7 
(a -1 = 2GeV). This gives us a great confidence in our approximate matching at one loop 
level, where such terms of wrong flavor channels, which enters at order g 2 , are neglected. 

From Figures |9| and [10], we also observe that the nearest neighboring data points have 
stronger correlation (less fluctuating) in the cubic wall source than in the even-odd wall 
source. The four times more floating point computation in the cubic wall source than in the 
even-odd wall source explains why the error bars of the data in Figure [1(] is about half of 
those in Figure Hence, we conclude that the even-odd wall source is equivalent to the 
cubic wall source from the standpoint of statistics. We believe that it is better to use the 
cubic wall source in weaker coupling simulations because the unwanted contamination from 
the degenerate and excited hadronic states becomes even more severe there. 



6 Data Interpretation 

Here, we interpret the numerical results in terms of the physics. First of all, we will address 
the technical questions about the improved wall source methods. Second, we will compare 
the numerical results of the one spin trace form and the two spin trace form and discuss 
the meaning of the consistency between the two formalism. Third, we will discuss the chiral 
behavior of Bk, and the chiral behavior of the individual components making up Bk (i.e. 
Ba, Bai, Ba2, By, By i, By 2 ) . In addition, we will discuss the effects of non-degenerate 
quark anti-quark pairs on the chiral behavior of these quantities. We will also discuss the 
chiral behavior from the standpoint of both 1/N C suppression and (partially quenched) chiral 
perturbation theory. Fourth, we will compare our numerical results with earlier works and 
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discuss the effect of quenching on Bk measurements. Finally, we will present our best value 
of Bk in the physical limit as well as our fitting procedure. 



6.1 Comparison of Wall Sources 

Numerical results for unrenormalized (naive, bare, or tree-level) Bk for various average quark 
mass, calculated both by the even-odd wall source and by the cubic wall source are shown 



in Figure [LI]. The values of Bk calculated by the two wall sources agree within errors. The 
error bars for Bk calculated by the cubic wall source are about half of those of the even-odd 
wall source. The cubic wall source method takes four times more computational time than 
the even-odd wall source method. We conclude that for Bk measurements at (3 = 5.7 the 
cubic wall source results are consistent with the even-odd wall source results. 



6.2 Comparison of One Spin Trace Form and Two Spin Trace 
Form 

The numerical results for unrenormalized Bk in both one spin trace and two spin trace 
methods are presented in Figure The numerical results for one-loop renormalized (tadpole 
improved through mean field theory 0, ||, ||, 0, ||]) B K both in the one spin trace 



and two spin trace methods are shown in Figure [13 
renormalized (RG improved (11 



The numerical results for one-loop 



34| ) Bk both in the one spin trace and two spin trace 
form are given in Figure |T3[ From Figures [12], [13] and |14], we notice that the results for the 
renormalized Bk calculated in the one spin trace method agree with those in the two spin 
trace method better than those for the unrenormalized Bk- 

Let us explain how we have obtained the renormalized coupling constant for the pertur- 
bative expansion. The one-loop renormalization of the AS* = 2 four-fermion operators on 
the lattice is explained in Ref. [11, [51, [5(J in detail. The detailed explanation of match- 
ing between the continuum and lattice observables at one loop level is also given in Ref. 
30|| . We discuss here how to obtain the coupling constant for the renormalization 



11, 31 



of the composite operators. For perturbative matching at one-loop level, one needs a well- 



defined coupling constant to use as the perturbative expansion parameter [RSI We have 



chosen the MS coupling constant at = ir/a scale as our perturbative expansion param- 
eter JHJ. There are two methods to obtain the MS coupling constant from the bare lattice 
coupling constant. One is a non-perturbative approach using tadpole improvement by mean 
field theory |33| and the other is a purely perturbative approach using renormalization group 



improvements |34| . We will refer to these two approaches to obtaining as the "tadpole 
improved" and "RG improved" methods. We use both methods to obtain the MS coupling 
constant at the renormalization scale it /a. 

The details of the non-perturbative approach of tadpole improvement by mean field 
theory is given in Ref 



IP 



The MS coupling constant (gf^) is related to 



the tadpole-improved coupling constant (g MF ) and the bare lattice coupling constant (g^) 
through mean field theory [O, EM: 



2 

9mf 



9o(a) 



ReUTiUn) 



(24) 
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9ms^ms) 



2 

9mf 



1 - Po9mf { 2 ln ( a VMs) + 2 ln 



-,9mf + 0{g MF ) 



'A 



Latt 



A 



A/5 . 



(25) 



where «7o(a) is the bare lattice coupling constant, U a is a unit gauge link plaquette, and 



0o 



11 

16^ 



2Nj 
33 



(26) 



For (3 = 5.7 and Nf = 2, we obtain the tadpole improved MS coupling constant: 

9mf — 1-8 2 

7T S 



2 

9ms 



a 



4^ - = 1.78 . 



(27) 



The coupling constant in Eq. (|27|) is used for tadpole improved matching between the 
continuum and the lattice observables. 

The details about a perturbative approach of RG improvement is explained in Ref. |34] . 
In this approach, the MS coupling constant (gjfg) is related to the bare lattice coupling 
constant (^(a)): 



9ms^ms) 



where t 



9o(a) 



t(3ogl(a) 



(28) 



M^ms) + ln 



'A 



Latt 



. A-MS , 



For j3 = 5.7 and Nf = 2, we obtain the MS coupling constant at /jtjjg = - scale: 

1.05263 

: 1.61 . (29) 



9ms I - 



The coupling constant in Eq. (^) is used for RG improved matching between the continuum 
and the lattice observables. 



6.3 Chiral Logarithms and Non-degenerate Quark Mass Pairs 

There have been a number of theoretical attempts to understand the chiral behavior of the 
K meson B parameters in terms of chiral perturbation theory p3| , |35| , [36| , |3~7| . In order to 
discuss the chiral behavior of the B parameters in an organized way, we need to consider a 
theory with four valence flavors: S and 5" both with mass m s as well as D and D', both 
with mass m^. Let the K° be the S*y 5 D pion and the K'° be the corresponding state with 
primed quarks (i.e. S'^^D' pion). Let us define B V1 , B V2 , Bai and B^ 2 as: 



B V1 = 4/r>,n i 577, Z rw l n\ /n i o. ~ r> I rsn\ ( 30 ) 



(K'° 


[S' allx D' h ][S bltl D a ] 


K°) 
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Bv2 

Bai 

Ba2 



(K'° 


[S' al ,D' a ][S bl ,D b ] 


K°) 






0)(0 




K°) 



(K'°\[S' alfll5 D' b }[S b ^ 5 D a ]\K°) 



(K 


s 




l lfil5 D> a \0)(0\S blfll5 D b \K°) 
[S> al » l5 D' a ][S bllll ,D b ]\K») 






0)(0 


S b i^D b 


K°) 



(31) 
(32) 
(33) 



where a and b are color indices. The Ba, By and B K are expressed in terms of Byi, B V2 , 
B A i and B A 2- 



B A 
By 
B K 



Bai + Ba2 
Byi + B V2 
By + Ba 



(34) 
(35) 
(36) 



Let us discuss the chiral behavior of the B parameters defined above {Bai, Ba2, Byi, 
By 2 ). In the vacuum saturation approximation with 1/Nq suppression, Bai — 0.25, Ba2 — 
0.75 and B V1 = B V2 = (obviously wrong!). For m u = m d ^ m s , chiral perturbation theory 
predicts the following results |I0, |23| : 



Byi 
By2 

Bai 

Ba2 



a 

(a_ - St) + a + A + (3i-^C + 

m k 

u 2 

(a+ - 5 2 ) + a_A + fc-^C + j 2 D 

mi 

P 2 

(a_ + Si) + a+A - f3i-^C - 

m\ 

n 2 

(a+ + 5 2 ) + ot-A - p 2 -^-jC - j 2 D 

mi 



(37) 
(38) 
(39) 
(40) 



where in the vacuum saturation approximation (tree-level chiral perturbation), a H 



~, 8i — \, S 2 = 1, fli = 1, (3 2 = h while 71 and 72 is not determined. A, C and D are 



defined as 
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m\ + m\ t 
2m\ 
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m K + m dd 

2m 2 K 



h{m dd ) 
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m\ + m 2 



2/i (m K ) - l2(m K ) 



2m\ 



'-h(m s 
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m 
A 2 " 



(44) 



m 2 / m 2N 



A 2 y (4tt/) 2 (4vr/) 2 V A ' 



(45) 



A is introduced as a momentum cut-off regularization for chiral perturbation theory. m 2 s = 
2/im s , m 2 dd = 2pm d and m 2 K = /i(m s + m d ) {i.e. 2m 2 K = m 2 s + m 2 dd ). Note there is 
no quadratic divergence in Byi, By2, Bai, and Ba2 even though ii(m) has a quadratic 
divergence. The quadratic divergences in C and D cancel out. The quadratic divergences in 
A is not a function of quark masses m s or m d . Hence, these quadratic divergences in A can 
be absorbed into the coefficients a_ and a + . As a summary of chiral perturbation theory, 
let us choose A = Airf and rewrite the results for B V1 , B V2 , B A1 , and B A2 as follows: 



B 



VI 



B 
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Bai 



B 
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where the constant terms proportional to are absorbed into <5j, and 



z = 



K 



(4tt/) 2 

ra s - rn d 
m s + m d 
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(50) 
(51) 



Here we keep full functional form with respect to e for theose terms of order ln(z) and z°. 
For those terms of order zln(z) and z, we do Taylor expansion with respect to e and keep 
only the terms of order e° and e 2 . The coefficients a iy rji, & and pi (i G {VI, V2, Al, A2}) are 
unknown. All the results in Eqn. ( |37] - |4Ci|) and Eqn. ([46]-||]) are calculated using full QCD as 
the fundamental theory. Note, By\, By2, Bai, and Ba2 have a branch point at e = ±1 which 
is non-singular. Let us set the domain ofeto— 1 < e < 1 and e G R. Byi, By 2 , Bai, and 
Ba2 have no singularity on this physical domain of e. Also note, Byi, By 2 , Bai, and Ba2 are 
even functions of e which means that although we switch m s with m^, the physics does not 
know it at all {i.e. it does not change). Hence, the power series expansion of B V1 , B V2 , Bai, 
and B A 2 with respect to e should have only even powers of e. The one-loop corrections to 
B V i, B V2 Bai and B42 include a term of the order hijnx) which is proportional to fa 
in Eqn. (|37| - |40|) . I 2 (rriK) /m 2 K is proportional to ln(m^) which is logarithmically divergent as 



ttik — > 0. These logarithmically divergent terms are called enhanced chiral logarithms [23 



and were originally noticed by Langacker and Pagels ||38|| . It is important to note that in full 
QCD, the enhanced chiral logarithms in Eqn. ( |37] - |40|) and Eqn. fl46| - |49|) are all proportional 

2 

to fa-^hiiTiK:), which depends only on the average mass of the quark antiquark pair. 

m K 

Let us adapt the above results to the (partially) quenched limit of QCD f23| . Here, 



quenched approximation means neglecting all the internal fermion loops and keeping only 
pure gauge interactions (i.e. sea quark mass is infinitely heavy), while the partially quenched 
approximation implies that the sea quark mass of the internal fermion loops is different 
from the valence quark mass. Both partially quenched and quenched approximations have 
additional infra-red problems which are absent in full QCD, since the sea quark fermion 
determinant can not regulate the infra-red pole singularity of the valence quark propagator. 
There are two important differences between (partially) quenched and full QCD. The first 
difference is that the meson eigenstates are not the same. The second difference comes from rj' 
loops, which can not contribute in full QCD simply because the 77' is too heavy. The diagrams 
of the hairpin type (77' loops) present in the (partially) quenched chiral perturbation for Bk 
vanish in the limit of m s = [p3|| . The logarithmically divergent terms which come from 
the hairpin diagrams are called quenched chiral logarithms [[24], 25, Here, the key point 
is that in (partially) quenched QCD, the quenched chiral logarithms in the B parameters, 
if present, must be a function of both average quark mass (m s + rrid) and mass difference 
(m s — ma) of the quark anti-quark pair, while the enhanced chiral logarithms common in 
both full and (partially) quenched QCD are not a function of quark mass difference but a 
function of only average quark mass. 

The quantity Bai i n Eqn. ( |39"D has a finite constant term which is a factor of around 1/3 
smaller than that in Ba2 in Eqn. (HO), while the enhanced logarithmic term in Bai is around 
3 times larger than that in -8,42- The enhanced logarithmic term in Byi is also around 3 times 
larger than that in By 2- Therefore, we have chosen Bai as a useful measurement adequate to 
observe both the enhanced chiral logarithms and, if present, the (partially) quenched chiral 
logarithms. 

We plot Bai and Ba2 with respect to the average quark masses in Figure |T3] and B V1 
and By2 with respect to the quark masses in Figure [16]. One can see a difference in the 
chiral behavior of Bai and By\ between the case with degenerate quark anti-quark pairs of 
mass: {0.01, 0.02, 0.03} and the situation with non-degenerate quark anti-quark pairs with 
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the mass pairs: {(0.004, 0.01), (0.004, 0.02), (0.004, 0.05), (0.01, 0.03)}. In order to more 
precisely see the effect of non-degenerate quark anti-quark pairs on B A i, we must interpolate 
in quark mass. Toward this end, we fit the data of the degenerate quark anti-quark pairs to 
the following function: 

B~ {mx) = Ai >o g (^) + M + ,o g (^) , ,52) 

and, second, subtract this form for the degenerate case from the non-degenerate data as 
follows: 



AB A1 (m K ) = ms + md (B A1 (m K , m s - m d ) - B A f enerate (m K )) , (53) 
m s -m d ^ 



where B A i(rriK, m>s — n^d) is the numerical data for the non-degenerate quark mass pairs. 
The covariance fitting results of the degenerate data on the jack-knifed ensembles are is 

B A T nerate (m K , a' 1 = 2.0GeV) = -0.683(47) log ( . m f + 1.305(97) 

+3Mis5) A ,o "(A) ^^^^ <*> 

B degenerate (m ^ ? Q -l = L g GeV) = _Q.300(25) log ( J^^) + 2.46(20) 

V( 47r .w / 

+6 .34(56)^lo g (^) (x » = 8.1X10-). (55) 



Using chiral perturbation theory, one can estimate the value of A\ in Eq. Q52]). From Eq. 

^ - -l^w^ • (56) 

where \i = m 2 K /{m s + ma). Using the results of our lattice QCD simulation (/i = 2.3a~ 1 ), 
this amplitude A\ is written in terms of a poorly known quantity (3\ ~ 1: A\ = — 1 1.6/?i . 
Using the value of meson mass and strange quark mass in the particle data book (m s = 
100 ~ 300MeV and m s jm d = 25), one can also obtain A\ = —(0.5 ~ 3.5)/?i in terms of a 
poorly known quantity j3i, where the range of values is chosen to reflect the large uncertainty 
in the chiral perturbation theory prediction. 

The numerical data for ABai{ttik) with respect to the average quark mass are plotted 
in Figure 0. From Figure [17], it is obvious that there is a strong sensitivity to the non- 
degenerate quark mass pairs (i.e. ABAi(mj() is not only a function of (m s + ma) but also a 
function of (m s — m d )). This additional divergence may be related to (partially) quenched 
chiral logarithms. It could also come from finite volume dependence on the lightest mass of 
the quark anti-quark pair. Something else might cause this unexpected divergence. The key 
point is that there is an additional divergence which is visible in our numerical simulations. 
At any rate, we would like to raise the following questions: 
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rf loop: Could the additional hairpin diagram in partially quenched chiral perturbation 
theory explain quantitatively the additional divergence which is sensitive to the non- 
degenerate quark masses? In fact, it is known that the contribution from rj' loops 



vanishes as m s — > 23, 10 



finite volume effect: Could the finite volume effect of Bai or By\ be sensitive to the 
non-degenerate quark masses? The small eigenvalues and their density in lattice QCD 
is regulated by the finite volume. Could these small eigenvalues be set up such that 
they are sensitive to the lightest mass of the non-degenerate quark anti-quark pair? 

fermion determinant in partially quenched QCD: In partially quenched QCD, the sea 
quark mass is different from the valence quark mass. In partially quenched QCD, the 
fermion determinant of the sea quark may suppress the coefficient of the logarithmic 
divergence more efficiently than the case of quenched QCD, which has no fermion 
determinants at all. Could one see the much larger effect of the non-degenerate quark 
anti-quark pair on Bai and Byi in quenched QCD than in partially quenched QCD? 



• scaling violation: The anomaly current in staggered fermion formulation x(7a*75 ® I)x 
is not a conserved current for finite lattice spacing. Hence, the corresponding pseudo- 
Goldstone pion x{l5 ® I)x nas a serious contamination of finite lattice spacing, which 
is supposed to vanish by flavor symmetry restoration in the continuum limit of a — > 0. 
This suggests that even in quenched QCD, 77' in staggered fermion formulation for 
finite lattice spacing may be much heavier than pseudo Goldstone pion x(ln ® £5)% 
and that its mass may have a scaling violation term of order a to make matters worse. 
How large is the scaling violation of rf mass as a function of quark mass? 

• other possibility: Is there something else related to the systematics of the non- 
degenerate quark masses on the lattice? 

One important thing is that this dependence on quark mass difference amplifies the en- 
chanced chiral logarithms rather than reducing them. We see this in Figure [IR O and 
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Ba and By are plotted in Figure [Tj| together with B K with respect to the average quark 
masses. One can notice that all these divergences related to both (partially) quenched chiral 
logarithms (if present) and enhanced chiral logarithms (present in Bai, Ba2, Byi and By 2) 
are canceled out in Bk which is finite in the chiral limit. 



6.4 Comparison with Earlier Work 

We now compare our numerical results of Bk with those of other groups as well as comparing 
our results at /3 — 5.7 (full QCD with Nf = 2) with the results at f3 = 6.0 (Quenched QCD). 
There have been two groups to calculate Bk at /3 — 5.7 with staggered fermions in full 
QCD. Kilcup |JT5 1 has calculated the unrenormalized (naive) Bk at /3 — 5.7 (16 3 x 32, 
full QCD with two dynamical flavors of a mass m sea a = 0.01, 0.015, 0.025). The lattice 
scale was a -1 = 1.9 ~ 2.0 GeV. The number of independent configurations was 50 and 
Bk measurements were done twice in different locations on the lattice for each individual 
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^valence 


Naive Bk 


this work (JN) 


this work (BS) 


Kilcup 


Ukawa et al. 


0.01 


0.600(27) 


0.595(22) 


0.658(18) 


0.69(2) 


0.02 


0.709(25) 


0.727(24) 


0.771(11) 


0.75(1) 


0.03 


0.768(23) 


0.791(21) 


0.818(09) 


0.79(1) 



Table 6: We compare our numerical results for naive Bk (cubic wall source, two spin trace 
form) with those of other groups (Kilcup and Ukawa et al). (3 = 5.7. m sea a = 0.01. JN 
implies that the errors are obtained through single-elimination jack-knife method. BS means 
that the data is analyzed by covariance fitting on the bootstrap ensembles. 



configuration, to make the total number of configurations equivalent to 100. Quark wall 
propagators were fixed in Landau gauge, and periodic boundary conditions in space and 
Dirichlet boundary conditions in the time direction were imposed. Ukawa et al. |fl9|, p0 



^■valence O 


Bai 


BA2 


this work 


Kilcup 


this work 


Kilcup 


0.01 


1.210(64) 


1.225(60) 


1.053(46) 


1.155(24) 


0.02 


0.565(23) 


0.575(15) 


0.860(30) 


0.932(13) 


0.03 


0.407(14) 


0.418(7) 


0.812(24) 


0.866(10) 



Table 7: We compare our numerical results of Bai and B^i (cubic wall source, two spin 
trace form) with those of Kilcup. (3 = 5.7. m sea a = 0.01. The errors are estimated through 
the standard jack-knife procedure. 



^valence O 


Byi 


Bv2 


this work 


Kilcup 


this work 


Kilcup 


0.01 


- 1.389(73) 


- 1.455(65) 


- 0.273(17) 


- 0.268(21) 


0.02 


- 0.638(26) 


- 0.664(21) 


- 0.0786(52) 


- 0.0713(53) 


0.03 


- 0.417(15) 


- 0.435(12) 


- 0.0342(22) 


- 0.0302(21) 



Table 8: We compare our numerical results of Byi and By 2 (cubic wall source, two spin 
trace form) with those of Kilcup. (3 = 5.7. m sea a = 0.01. The errors are estimated through 
the standard jack-knife procedure. 



studied Bk at (3 = 5.7 on a lattice of size 20 3 x 20 (duplicated in the time direction) with 
two flavors of dynamical staggered quarks of mass m sea a = 0.01. The lattice scale was 
a = 0.085 ~ 0.09 fm (a -1 = 2.2 ~ 2.4 GeV). Both Landau gauge operators and gauge- 
invariant operators were used. Quark propagators were calculated with Dirichlet (periodic) 
boundary condition in the time (space) direction. 

The differences between our numerical simulation and those of other groups' are the 
lattice size, the boundary conditions on the quark propagators in the time direction and 
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^valence " 


Naive Bk 


full QCD (N f = 2) 


quenched QCD 


this work (JN) 


this work (BS) 


Sharpe et al. 


Kilcup 


Ukawa et al. 


0.01 


0.600(27) 


0.595(22) 


0.68(2) 


0.697(29) 


0.69(2) 


0.02 


0.709(25) 


0.727(24) 


0.73(1) 


0.749(16) 


0.74(1) 


0.03 


0.768(23) 


0.791(21) 


0.78(1) 


0.777(13) 


0.78(1) 



Table 9: The full QCD calculation of naive (unrenormalized) Bk (cubic wall source, two 
spin trace form) is compared with that of quenched QCD (Sharpe et a/., Kilcup and the 
Ukawa et al. group). JN means that the errors are estimated through the standard jack- 
knife procedure. BS means that the data is analyzed by covariance fitting on the bootstrap 
ensembles. 



the color summation over the meson wall sources. We are not summing the three values of 
the color index for the meson wall sources on the individual configuration sample. Instead, 
we choose a different color index for the meson wall sources in each measurement (in other 
words, color indices of meson wall sources are spread over the configuration samples with 
equal statistical weight instead of being summed on each configuration sample). The Bk 
results of this work, Kilcup and Ukawa et al. are summarized in Table ^| The Bai, Ba2, Byi 
and Bv2 of this work and Kilcup are compared in Tables ^ and |8|. We believe that the 
differences in Tables || |7] and |8] can be reasonably explained by the different lattice size, the 
boundary conditions on the quark propagators, color summation of the meson wall source, 
poor statistics, the uncertainties in the lattice spacing, etc.. Hence, we conclude that all the 
measurements in Tables ||, [7] and [8] are consistent with one another. 

It is also important to compare the results of full QCD (Nf = 2) with those of quenched 
QCD. There have been three groups to calculate Bk at /3 — 6.0 (a" 1 = 2.0 ~ 2.1 GeV) 
in quenched QCD. This corresponds to (3 = 5.7 in full QCD. The results of this work (full 
QCD, = 5.7, 16 3 x 40), Sharpe et al. [pj, |I§ (quenched QCD, (3 = 6.0, 24 3 x 40), Kilcup 
P, |T| (quenched QCD, (3 = 6.0, 16 3 x 40) and Ukawa et al. p], |2| (quenched QCD, 
(3 = 6.0, 24 3 x 40) are compared in Table |9|. From Table |^, we conclude that the dependence 
of Bk on m sea is too weak to detect in our numerical simulations (in other words, the effect 
of quenched approximation is less than 15 % in Bk measurements). This has been predicted 



by the quenched chiral perturbation theory in the limit of m s — > [10, 



6.5 Fitting Procedures for Bk 

Bk describes K° — K° mixing at the energy scale of about 500 MeV (i.e. in the low energy 
limit of QCD dynamics). It is not known how to calculate the dependence of Bk on the 
valence quark mass directly from the QCD Lagrangian. For this reason we adopt the chiral 
effective Lagrangian (equivalent to current algebra), which is valid in the energy region 
below the p meson mass. This chiral effective lagrangian is not a cure-all solution to the low 
energy dynamics of QCD. However, it gives us a reasonable guide to understand the leading 
chiral behavior of QCD. We will now use the predictions of the chiral effective Lagrangian 
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to interpolate between our B K results to make a B K prediction for physical quark masses. 

The corrections from chiral perturbation theory to M.k in full QCD were calculated in 
Ref. p3| , |35| , p6| , 57]. The results of that calculation |23|] [m u = m d ^ m s ) in full QCD were 
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where e is defined as 



e = 



(57) 



(58) 



Ci, C2 are unknown coefficients, but can be determined by numerical simulation on the lattice. 
The results from quenched chiral perturbation theory |23| [m u = m d ^ m s ) are 
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where quenched chiral perturbation |23[ predicts 



5 = 



Ami 



0.2 



(59) 



(60) 



c[ and c 2 are unknown coefficients. The term proportional to 5 is a contribution of rj' loops 
appearing in quenched QCD, which are absent in full QCD. Because of this term, j^^ enched 
has a singular branch point at e = ±1. Note, B^ 11 does not have any singular branch 
point. Equations fl5"T|) and fl5"9|) are the theoretical predictions for B K as a function of the 
light quark masses. 

Let us choose our fitting function on the basis of the predictions of (quenched) chiral 
perturbation theory in Eq. (|57|) and (|59|). Parts of our numerical simulation of Bx are 
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classified in the category of partially quenched QCD (Nf = 2 but m sea ^ m valcnce ). We 
notice that the coefficients of those terms proportional to e 2 m? K ln(m^-) in Eqn. (|57|) and 
Eqn. ( |59| ) are different from each other, which implies that the coefficient of these terms 
should be determined by our numerical data. 

We have tried a linear fitting function (BK(m q ) — a.\ + ci2m q a) even though the linear 
term is a next to leading order correction from the chiral perturbation (as can be seen in Eqn. 
(|5"7|) and QoTj|)). The x 2 /(d.o.f.) for the linear covariance fitting on the jack-knifed ensembles 
of Bk data calculated in the two spin trace form with the cubic wall source method is 
17.2(37). This implies that the fitting is poor and that we need an additional term (e.g. 
m q a\n.(m q a)) to fit the data. 

Hence, we choose the first fitting function as 

-Bi^(m g ) = «i + a2m q a\n(m q a) + a^,m q a . (61) 

where (i — 1, 2, 3) are unknown coefficients to be determined and m q a = \(m s + m^). In 
this first fitting function, we neglected the effect of non-degenerate quark mass pairs (i.e. 
terms proportional to e 2 are neglected). From the chiral perturbation theory Eqn. ([57D and 



(|59|), the predicted value of the ratio a>2 /a% is 



- = ( 3 + A )77^' (62) 

where < A < 1. Using our lattice QCD simulation results for //(// = 2.3a -1 ), the ratio 
02/ ci\ is written in terms of poorly known quantity A: 

— = 40 + 13A . (63) 



s — 



Using the value of K meson mass and strange quark mass in the particle data book (m 
100 ~ 300MeU and m s /md = 25), one can also obtain the ratio a 2 /a\ 

— = (7.0 ~ 23) + (2.3 ~ 7.6)A, (64) 

in terms of a poorly known quantity A, where the range of values is chosen to reflect the 
large uncertainty in the chiral perturbation theory prediction. The Bx results of the above 
3 parameter covariance fitting on the bootstrap ensembles are summarized in Table [10] (un- 
renormalized) and |TI] (tadpole-improved renormalized) . From Table [Tl] and O, we observe 
that the renormalization with tadpole improvements makes BK(m q ) in agreement between 
the one spin trace form and the two spin trace form once they are obtained with the same 
source method and covariance-fitted in the same range. The covariance fitting results for 
the ratio 0:2/01 were not consistent among the various measurements, mainly because the 
3 parameters have a wide domain to fit the 7 B^ data points as a function of the average 
quark mass. This gave us a motivation for the second fitting trial function, which will have 
only two parameters. 

Our second fitting trial function is chosen such that the ratio Oi/a^ is fixed to the case 
of the degenerate quark anti-quark pair (e = 0). The reason is that this ratio is universal for 
the degenerate quark pair regardless of quenched approximation and that our e values are 
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small enough to take into consideration later as a perturbative expansion parameter. The 
second fitting trial function is 



B K {m K ) = q;i 
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(65) 



where aj (i = 1,2) are unknown coefficients. We take the effective mass (mx) of fx ((0 | 
| K )) measurements in Table f| as a definition of mx in Eq. (|65|) . In this second fitting 
trial function, we set the coefficient of the leading term in the chiral perturbation expansion 
to the theoretically expected value in Eq. ([57], [5^) and also we again neglect the effect of 
non-degenerate quark mass pairs. When we set a~ l = 1.8 GeV from the p meson mass and 
choose f n = 93 MeV, the covariance fitting results for the two spin trace form (cubic wall 
source) are summarized in Table ^ and |I3[ From Table [T2| and |T^, we notice that the fitting 
results for the tadpole- improved renormalized Bx in various types of the measurements are 
in good agreement with one another, while the fitting results for the unrenormalized Bx are 
not consistent between the one spin trace form and the two spin trace form. 

Let us discuss how we can detect the effect of non-degenerate quark anti-quark pairs 
on Bx (i-e. the dependence of Bx on e). The strategy is the following. First, we divide 
the numerical data in two parts: one part (we call this the degenerate part) contains only 
the numerical data for degenerate quark mass pairs {(0.01, 0.01), (0.02, 0.02), (0.03, 0.03)} 
and the other part (the non- degenerate part) contains only the data for non-degenerate 
quark mass pairs {(0.004, 0.01), (0.004, 0.02), (0.004, 0.05), (0.01, 0.03)}. Next, we fit the 
degenerate part to the second trial function in Eqn. ([65]) which are supposed to be exact for 
the degenerate quark mass pairs up to the given order in the chiral perturbative expansion. 
For the degenerate part, 



B 



degenerate 
K 1 



m K 



1-3 



in 



K 



In 



m 



K 



+ a 2 



in 



K 



(47TA) 



(66) 



Now we define a function which represents the effect of non-degenerate quark mass pairs: 



AB K {m K ) = 



B K (m K ,e) 



^degenerate / \ 

B K {m K ) 



(67) 



where Bx{mx,£) are our numerical data with non-vanishing e 2 = {m s — md) 2 /(m s + m^) 2 
and i?J gcneratc (mx) represents the fitting function for the degenerate part, given in Eq. (|66|) . 
Finally, we try to find a functional form, if possible, to fit ABx{mx) numerical data to. In 
Figure ^ (two spin trace form, cubic wall source), Figure ^ (two spin trace form, even-odd 
wall source) and Figure ^l| (one spin trace form, even-odd wall source), we plot ABx with 
respect to the average quark mass, only for the non-degenerate quark anti-quark pairs. As 
you can see in the figures, it is hard to find a functional form which can explain all of the 
data. From these figures, we notice that the dependence of Bx on e 2 is extremely small 
(especially in the domain near the physical K meson mass). Hence, we conclude that we 
could not detect any significant effect of non-degenerate quark mass pairs on Bx within the 
precision of our numerical study. 

As a conclusion to this section, let us present our best value of Bx- Since the data of the 
cubic wall source method has better statistics than that of the even-odd wall source method, 
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Trace Form (Source) 




a 2 


a 3 


X 2 /d.o.f. 


B K {m q ) 


2TR (Cubic) 


0.336(65) 


-11.2(33) 


-24.1(93) 


1.05 


6.600(21) 


2TR (Even-Odd) 


0.346(89) 


-10.1(39) 


-21.9(112) 


1.87 


6.375(29) 


1TR (Even-Odd) 


0.477(119) 


-9.35(482) 


-21.9(137) 


0.47 


7.270(30) 



Table 10: Covariance fitting of unrenormalized B K with 3 parameters: Here 2TR and 1TR 
represent the two spin trace form and the one spin trace from respectively. Cubic and Even- 
Odd imply the cubic wall source and the even-odd wall source respectively and B K {m q ) 
means interpolation to the physical quark mass. 



Trace Form (Source) 


a\ 


«2 


«3 


X 2 /d.ol. 


B K (m q ) 


2TR (Cubic) 


0.372(62) 


-9.83(32) 


-21.6(95) 


1.06 


0.655(21) 


2TR (Even-Odd) 


0.416(86) 


-7.43(37) 


-15.9(102) 


2.03 


0.636(28) 


1TR (Even-Odd) 


0.408(112) 


-8.30(480) 


-19.2(122) 


0.52 


0.633(27) 



Table 11: Covariance fitting of tadpole-improved renormalized Bk with 3 parameters at 
the scale of fi = n/a: Here 2TR and 1TR represent the two spin trace form and the one spin 
trace from respectively. Cubic and Even-Odd imply the cubic wall source and the even-odd 
wall source respectively and BKijn q ) means interpolation to the physical quark mass. 



Trace Form (Source) 




« 2 


X 2 /d.o.f. 


B K (m K ) 


2TR (Cubic) 


0.293(14) 


0.397(66) 


0.897 


0.638(21) 


2TR (Even-Odd) 


0.289(21) 


0.334(82) 


1.53 


0.619(28) 


1TR (Even-Odd) 


0.346(24) 


0.191(84) 


0.70 


0.702(35) 



Table 12: Covariance fitting of unrenormalized Bk with 2 parameters: Here 2TR and 1TR 
represent the two spin trace form and the one spin trace from respectively. Cubic and Even- 
Odd imply the cubic wall source and the even-odd wall source respectively and B K {m K ) 
means interpolation to the physical K meson mass. 



Trace Form (Source) 




a 2 


X 2 /d.o.f. 


B K (m K ) 


2TR (Cubic) 


0.301(14) 


0.300(67) 


0.89 


0.636(21) 


2TR (Even-Odd) 


0.297(21) 


0.255(83) 


1.83 


0.619(28) 


1TR (Even-Odd) 


0.300(21) 


0.184(81) 


0.47 


0.612(29) 



Table 13: Covariance fitting of tadpole-improved renormalized Bk with 2 parameters at 
the scale of /i = ir/a: Here 2TR and 1TR represent the two spin trace form and the one spin 
trace from respectively. Cubic and Even-Odd imply the cubic wall source and the even-odd 
wall source respectively and B K {jn K ) means interpolation to the physical K meson mass. 
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we shall quote the covariance fitting results of the cubic wall source data as our best value. 
For the second trial function, the uncertainty in the lattice scale a and the decay constant 
fir (fx) produces uncertainty in Bk{ttik) of the same order of magnitude as the statistical 
uncertainty, while covariance fitting to the first trial function is not as much sensitive to the 
lattice scale a and the decay constant f w . Hence, we choose the results of the covariance 
fitting to the first fitting trial function with three parameters on the bootstrap ensembles as 
our best value. Interpolated from the first fitting trial function, the physical results are 

unrenormalized B K (m K ) = 0.660(21) 
renormalized (N.D.R.) B K (m K ,fi = -) = 0.655(21) , (68) 

(X 

where N.D.R. implies naive dimensional regularization scheme and the errors represent purely 
statistical uncertainty. Here the physical results imply the Bk values for the physical kaon 
mass (rriK = 497.7 MeV). Here, we have completely neglected the systematic errors related 
to the scale (a) uncertainty, the coupling (flfjg) uncertainty, the contamination from the 
unwanted operator mixing, and the contamination from unwanted hadronic states which 
can couple to the operators used. We also could not control finite volume effects, finite 
temperature effects, or finite lattice spacing effects. 

7 Conclusion 

Here, we summarize what we have learned through the numerical simulation of Bk and what 
needs further investigation in the future. 

The results for B K from the improved wall source (cubic wall source) were in good 
agreement with those of the conventional even-odd wall source. It is shown that the cubic 
wall source suppresses the contamination from the wrong flavor channels efficiently However, 
the cubic wall source takes four times more computational time than the even-odd wall 
source. In the limit of a — > 0, the SU(4) flavor symmetry is supposed to be recovered and 
so there will be a serious contamination from pseudo Goldstone pions with various flavor 
structures and the various p mesons. The cubic wall source is quite promising in the weak 
coupling region to exclude the contamination from unwanted hadronic eigenstates. 

We can transcribe the continuum AS* = 2 operator to the lattice with staggered fermions 
in two different ways. Theoretically, both formalisms of operator transcription must be 
equivalent to each other in the limit a — > 0. The numerical results in the one spin trace 
formalism were consistent with those in the two spin trace formalism after the proper renor- 
malization (with either tadpole or RG (MS) improvement). We have learned how important 
the proper renormalization is, as well as the careful choice of the coupling constant for the 
perturbative expansion. 

We have tried to understand the effects of the non-degenerate quark masses on B K and 
the individual components making up Bk- The effects of the non-degenerate quark anti- 
quark pairs on Bk were too small to observe within the precision of our numerical simulation 
(especially near the region of physical K meson mass) . Why this effect is so small needs more 
careful theoretical investigation. Chiral perturbation theory suggests that B^i (axial part 
of B K with one color loop) is the best observable to detect the enhanced chiral logarithms 
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which are not expected to be a function of quark mass difference. We observed an additional 
divergence which depends on the quark mass difference. This additional divergence needs 
more thorough investigation and understanding. We wonder whether partially quenched chi- 
ral perturbation can explain this additional divergence, or how much finite volume effects on 
-B41 depend on the lightest mass of the non-degenerate quark anti-quark pair. Qualitatively, 
chiral perturbation theory is consistent with our numerical work. However, more theoretical 
research on the (partially) quenched chiral perturbation and its quenched chiral logarithms is 
necessary to see whether the (partially) quenched chiral perturbation can explain the effects 
of the non-degenerate quark antiquark pairs on B^ and its individual components. 

We could not observe any dynamical fermion effect on B K . It is difficult to understand 
why internal fermion loops are so unimportant for B K , since the Dirac eigenvalue spectrum 
of quarks is supposed to be completely different between quenched QCD and full QCD. This 
also needs further theoretical understanding. 

Through the numerical study of Bk in this paper, we have learned that the cubic source 
method is promising for the weak coupling simulation and that the one spin trace formalism 
is consistent with the two spin trace formalism. It is true that lattice QCD results for B K 
are more solid and believable after this work. However, it is also true that there are many 
details which need more thorough investigation and understanding. 
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Figure 1: Vacuum saturation with respect to time, maa = m s a = 0.01. Calculated with the 
even-odd source method in the two spin trace formalism. 
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Figure 2: Vacuum saturation with respect to time, rrida = m s a = 0.01. Calculated with the 
cubic wall source method in the two spin trace formalism. 
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Figure 3: Mk with respect to time, m^a = m s a = 0.01. Calculated with the even-odd wall 
source method in the two spin trace formalism. 
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Figure 4: Mk with respect to time, m&a = m s a = 0.01. Calculated with the cubic wall 
source method in the two spin trace formalism. 
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Figure 6: Unrenormalized Bk with respect to time, m^a = m s a = 0.01. Calculated with 
the even-odd wall source method in the two spin trace formalism. 
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Figure 7: Unrenormalized Bk with respect to time, mja = m s a = 0.01. Calculated with 
the cubic wall source method in the two spin trace formalism. 
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Figure 8: Unrenormalized Bk with respect to time, mja = m s a = 0.01. Calculated with 
the even-odd wall source method in the one spin trace formalism. 



38 



0.04 



CD 

Ctf 
O 

o 
> 

H — 

O 

% 
CD 



CD 




15 20 

Time 



Figure 9: Unrenormalized B K with the wrong flavor structure (((V+ A) ®S) 2TR with respect 
to time, mda = m s a = 0.02. Calculated with the even-odd wall source method in the two 
spin trace formalism. 
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Figure 10: Unrenormalized B K with the wrong flavor structure (((V + A) <g> S) 2TK with 
respect to time, m^a = m s a = 0.02. Calculated with the cubic wall source method in the 
two spin trace formalism. 
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Figure 11: Unrenormalized with respect to average quark mass. The filled circles rep- 
resent the data from the cubic wall source. The empty squares represent the data from the 
even-odd wall source. All the data are obtained using the two spin trace formalism. 
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Figure 12: Unrenormalized B K with respect to average quark mass. The filled circles repre- 
sent the data from the one spin trace form. The empty squares represent the data from the 
two spin trace form. Both data are obtained using an even-odd wall source. 
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Figure 13: Tadpole-improved renormalized Bk with respect to average quark mass, fi = ir/a. 
The filled circles represent the results from the one spin trace form. The empty squares 
represent the results from the two spin trace form. Both results are obtained using an 
even-odd wall source. 
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Figure 14: RG improved (MS) renormalized Bk with respect to average quark mass, fi = 
it /a. The filled circles represent the results from the one spin trace form. The empty squares 
represent the results from the two spin trace form. Both results are obtained using an 
even-odd wall source. 
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Figure 15: Bai, Ba2 with respect to average quark mass. Three data points with average 
quark mass G {0.01,0.02,0.03} correspond to the degenerate quark masses. The other four 
data points correspond to the non-degenerate masses. The data are obtained using the two 
spin trace formalism with the cubic wall source method. 
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Figure 16: Byi, By2 with respect to average quark mass. Three data points with average 
quark mass G {0.01,0.02,0.03} correspond to the degenerate quark masses. The other four 
data points correspond to the non-degenerate masses. The data are obtained using the two 
spin trace formalism with the cubic wall source method. 
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Figure 17: ABai with respect to average quark mass. The quark mass pairs are (0.004, 
0.01), (0.004, 0.02), (0.01, 0.03) and (0.004, 0.05). The data are obtained using the two spin 
trace formalism with the cubic wall source method. 
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Figure 18: Unrenormalized B v , B K and Ba with respect to average quark mass. The data 
are obtained using the two spin trace formalism with the cubic wall source method. 
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Figure 19: AB K with respect to average quark mass. The data are obtained using the two 
spin trace formalism with the cubic wall source method. 
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Figure 20: AB K with respect to average quark mass. The data are obtained using the two 
spin trace formalism with the even-odd wall source method. 
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Figure 21: AB K with respect to average quark mass. The data are obtained using the one 
spin trace formalism with the even-odd wall source method. 
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